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ABSTRACT 

In this paper , a “Proportional-Integral-Derivative” (PID) controller used in a aircraft pitch control system is 
optimally tuned to control the output response so as to obtain a much more stable and optimized transient response is 
studied. The objective is to obtain a cohesive, sturdy and controlled system by tuning the PID controller using PSO 
algorithm. The obtained result is compared against the results obtained by re-tuning the same system with standard PID 
tuning algorithm and proves to give more suitable results. Results obtained establishes that , if we tune the PID controller 
of the given system using PSO algorithm we get reduced percentage of overshoot reduced rise time, settling time, peak 
and hence system is more stable and less stagnant with PSO than the same system tuned with standard PID tuning 
algorithm. 
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INTRODUCTION 

Several control systems of aircraft and other control industries use scholastic controller like PID to upgrade 
the systems features and dynamic conduct and to improvise the stability assay and performance dynamics of the 
aircraft. The flight envelope safety is being enhanced by optimal tuning of parameters of PID controller, for 
aircraft’s dynamics of pitch control. The wide significance of the PID controllers can be attributed, to their 
admirable performance, for a wide operating condition ranges, and simplicity of functioning, which lets engineers 
operate them in a simple, straightforward manner and commonness, with which it is anticipated amongst researchers 
and professionals in the area of process control industries. PID controller calculates the net value of inaccuracy, 
which is the discrepancy in output result of the system and optimum reference value. It reduces the error value, by 
modulating the inputs of the pitch control system. 

PID controller gains are 3-term control parameters - the Proportional, Integral and Derivative parameters 
denoted as K P , IQ, and K d . The output of the controller is presented in terms of error in value of the controller, the 
extent to which it deviates from the standard result achieved, calculating of rate of oscillation of the system. 
The structure of PID controller and is represented by- 

X (s) = Kp+Krfl/s) +K ds = K P [1+ (1/Ts)+T ds ] (1) 

Where K P is proportional gain parameter, IQ is integral parameter gain, K d derivative gain parameter, T r is 
time constant integral and T d is time constant derivative. A proportional controller (K P ) reduces the rise time and it 
also reduces, but doesn’t get rid of the steady-state error. An integral parameter (IQ) eliminates the steady-state 
error, for a constant or step input, but it may slow down the transient response for the system. A derivative control 
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(K d ) increases the stability of the system, reduces the overshoot percentage, and improves the transient response”. 
These correlations are not precisely accurate, because K P , Kj& Revalues, depend on each other and modifying one of these 
variables, can change the values of others. 

Despite of widespread application of PID controller’s one specific limitations is that there is no efficient method 
to tune this type of controller. Various methods have been proposed over the years to tune the PID controllers. Some of 
them will be discussed later on in this paper. 

As mentioned before for this paper the swarm optimisation algorithm of PSO is used to tune PID gains/parameters 
(Kp, Ki, Kd), which in turn is used to optimize the longitudinal pitch control parameters, for an aircraft pitch control 
system. The aim of this paper is to establish that, if we tune the PID controller using PSO technique, we get less overshoot, 
less rise time, less settling time, less peak time and hence, system is less stagnant and much more stable, compared to Mat 
lab auto tuning algorithm. 

Pitch Control System Using PID 

Over the years, flight control systems are being engineered, by employing mathematical specimen of the airplane 
linearized for several dynamic flight parameters modified, with the operating situations. The purpose of the work 
developed is, to modulate the pitch rate of an airplane for pitch control, in order to give the system better stability related 
characteristics, when the airplane is upward nosed or downward nosed. The system for control of pitch rate discussed in 
this paper, is demonstrated in Figure. 1. Here, “X b , Y b and Z b ” denote the aerodynamics force elements. “0, O and 8 e ” 
denote the orientation of aircraft, w. r. t the earths-axis model and the angle of deviation of the elevator. 


Figure 1: Aircrafts dynamic Pitch control Figure2: velocity, force and moment components 

Model Parameters [7] in rigid body coordinate axes [7] 

Figure.2, shows the velocity, force and moment constituents in the body fixed coordinate of aircraft system. The 
aerodynamics moment constituents for pitching, yawing and rolling axes are denoted as M, N and L The term, q, r and p 
represent the angular rates with respect to pitching, yawing and rolling axis while term u, v, w denote the velocity 
components of rolling, pitching and yawing axes. Angles a and P denote the angle of attack and sides lip values 
respectively 

The assumptions considered for the given aircraft model used in this paper are 
The assumptions and equations considered for the modelling process are 

• A constant state of cruising at steady altitude and velocity, therefore the aircrafts drag and thrust value nullify 
each other and lift generated and weight balances one another 

• pitch modification doesn’t affect velocity of the airplane in any manner 
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• Rolling rate P ^ 6 

• Rate of yaw, ^ — 8 

• Pitch Rate r #C& - 0S& 

• Pitch Angle, @ <t> 1 A <i> 

• Angle of Roll & ~ P + c lS <t> T 0 + rC 0 T e 

• Angle of Yaw O ) sec & 

We refer to the Figure 1 and Figure 2, from which dynamic equations of force and moment are determined. 

X - mgS d = m(u + qv- rv) 

Z + mgC 0 C tp = m(w+ pv-qu) 

M = 1 y q + rq(! x - /.) + l^p 2 - r 2 ) (3) 

Equations of [1] [2] and [3] are being linearized with “small disturbance theory”. The equations are replaced by a 
particular standard variable, with a perturbation or disruption, as it is shown below. 

u = u Q 4- Aw v = v Q + Av w= w Q + Aw 

P = Po +A P q = r = r 0 +Ar 

X = X a +AX M -M 0 +MY Z = Z 0 + AZ 
3 — $ D + A $ 


For our convenience, we assume the standard conditions of flight to be balanced and the propulsive forces as 
constant. This implies that, " v 0 = po = qo = r 0 = O 0 =\j/ 0 = w 0 = 0". From linearization the eons [4], [5] and [6] are being 
given below 



An - X w Aw+(g cos 0 o )A 6 = XS e AS e 


( 4 ) 



( 5 ) 

( 6 ) 


We manipulate values of eqns [4], [5] and [6] and substitute the value of parameters of the derivatives of 
longitudinal stability, thus we obtain this particular transfer function, for the rate of change of pitch, to the change in angle 
of deflection of the elevator is being denoted as [7] 
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W = ~(M & +M a Z & lu,)s-(MA /»„ -M & ZJu 0 ) 

M( s ) r ~(M q +A/ 4 +ZJu li )s+{Z a M q I w 0 -M a ) (?) 

The transfer function for the rate of change of angle of pitch, to the rate of deflection of the elevator can be 
acquired, from the pitch rate variance to the change in angle of elevator in the following way. 

Aq = Ad 

( 8 ) 

Aq(s) = sA0(s) (9) 

A0(s) _ 1 Aq(s) 

A S c (s) s A 9(s) (1Q) 

Thus the transfer function for the given airplane pitch control system is achieved in eqns [11] and [12], 
respectively. 


A 0(s) _ 1 -(M & + M dr Z & /u 0 )s-(M a Z & /u 0 -M & Z a /u 0 ) 
A4(s) 5 * s 2 -(M q +M a + Z a /u 0 )s + (Z a M q /u n - M a ) 


A 0(s) 11,73045 + 22.578 

AS e (s) ~ s 3 +4.9616s 2 + 12.941s 

( 12 ) 

This transfer function block, along with the actuator is implemented using PID to give the final model, as it is 
shown in the below figure 



Figure 3: Pitch Control System Using PID Controller [6] 

A proportionality parameter gives the signal error value, through the constant gain aspect indicators. An integral 
parameter helps decreasing steady-state error and a derivative parameter improves the transient feedback of the aircraft 
system. A collective PID controller generally gives more desirable solution than independent of P, I and D parameters. 
Choosing of gain value, for the PID controllers can be obtained from a separate closed loop method of tuning 


• An aircraft pitch control system provides control action for longitudinal pitch modulation of the aircraft. 


• The control is provided in the form of pitch rate and pitch angle at a specific point of flight 


PID Tuning Methods 
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The various primary tuning methods employed these days are “Ziegler Nichols” method, “Modified Ziegler 
Nichols” method, “Asstrom and Haggland” method, “Tyrus Leuben” methods etc. 

PID Tuning Using PSO for Pitch Control System 

PSO is swarm intelligence based indeterminate optimization method constructed and influenced by social 
behaviour of flock of bird or school of fish. PSO simulates the behaviours of these social systems. For PSO, we consider a 
particular solution as a "bird" in the search territory. We name it "particle". Every single particle has a particular value of 
fitness which is calculated by the fitness function which we want to optimize; every single particle has a particular velocity 
which directs its flying. The particles fly through the problem space by following the particle with the best fitness function 
PSO computes a group of particles (solutions) in random space and after that inspects for optimum value by revising he 
values through multiple iterations. For each iteration, every single particle is updated by pursuing the two "best" results. 
The 1 st value is the best result (fitness), it has attained so far. (That value of fitness is saved.) This value is known as 
“pbest”. Another value "best", that is pursued by the particle swarm optimizer is the best result, being attained so far by 
any particle in the whole group. This value is a global best and named “gbest”. If one particular particle takes portion of the 
population as its metaphoric acquaintances, the optimum value is a local best and is called 1 best. 

After the two best values are obtained, the velocity and position of each particle is updated with the given 
equation (a) and (b). 

• V P - V P + A1 * rand() * (pbest - current value) + A2 * rand() * (gbest-current value)(i) 

• Current value = current value + V p (ii.) 

Vp is the velocity of the particle; current value is the present particle fitness (solution), “pbest” and “gbest” are 
defined as given earlier, rand () is a random value between (0,1). Al, A2 are learning factors which in several cases are 
Al= A2 = 2 



Figure 4: Pitch Control System with PID Controller using PSO [4] 

This is the final model developed for the system using PSO. As shown above here PSO algorithm is used as a 
function block with a function caller to tune the PID controller and thus give a tuned value of Kp, Ki and Kd, which is then 
fed into the system to give us the desired transient response. 

RESULT & DISCUSSIONS 



Figure 5: PSO Tuned Response Vs Matlab Auto Tuned Response 
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The following graph displays the final transient response and compares obtained PSO response with Matlabauto 
tuned response 

From the above Amplitude vs. Time graph we can observe that the PSO tuned response for the aircraft pitch 
control looks way more stable than the Matlabauto tuned response, for the same system. This, in particular, is an extremely 
important result because, a more stable system would mean enhancement in the manoeuvrability or control of the aircraft, 
in control of pitch rate. Having a good manoeuvrability is an extremely vital part for an aircraft, during its flight because; 
the safety of the flight directly depends on it. Moreover, the efficiency of the aircraft also increases due to a more stable 
pitch control system, thus reducing the cost of operations. 

From the two transient responses above, we obtain the performance parameters - Rise time, Peak time, Settling 
time, Maximum Overshoot percentage. We write down these parameters in Tabular form for analysis and then compare 
them. 

PID Controller Parameters 


PID Parameters 

Auto tuned 

PSO tuned 

K P 

7.655 

0.63 

K, 

11.82 

0.0504 

K d 

1.189 

1.9688 


Performance Parameters 



Auto Tuned 

PSO tuned 

Rise Time 

0.0885 

0.0504 

Settling Time 

1.15 

0.329 

Overshoot % 

9.52 

4.31 

Peak 

1.1 

0.875 


CONCLUSIONS 

We can see from the graph and tabular results above that the transient response for the pitch control system model 
with PSO gives much more stable and optimized response compared to the auto tuned response obtained for the same 
system as all the performance parameters obtained-Rise time, Settling time, Overshoot percentage and peak are much less 
for PSO tuned response than Mat lab auto tuned response. 

The proves that, for the given aircraft pitch control the results obtained from PSO algorithm are much more 
efficient and adequate than those obtained from Matlabauto tuning algorithm. 
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